Detecting monopoles on the lattice 
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We address the issue why the number and the location of magnetic monopoles detected on lattice 
configurations are gauge dependent, in contrast with the physical expectation that monopoles have 
a gauge-invariant status. By use of the non-Abelian Bianchi identities we show that monopoles are 
gauge-invariant, but the efficiency of the technique usually adopted to detect them depends on the 
choice of the gauge in a well understood way. In particular we have studied a class of gauges which 
interpolate between the Maximal Abelian gauge, where all monopoles are observed, and the Landau 
gauge, where all monopoles escape detection. 
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I. INTRODUCTION 

Monopoles may play a fundamental role in gauge theo- 
ries, their condensation in the vacuum being responsible 
for dual superconductivity and color confinement QUI]- 

The correct way to detect monopole condensation is 
to measure the vacuum expectation value of an operator 
fj, carrying nonzero magnetic charge: if (/x) =/= vacuum 
is a dual superconductor [4j. In the deconfined phase 
(//) =0. Much work has been done on this line [a-lllj. 

Much more activity, however, has been devoted to the 
observation of monopoles in lattice configurations, cither 
to look for an effective monopole-action, based on the 
empirical observation that monopoles could be the dom- 
inant degrees of freedom of the theory (monopole domi- 
nance) [L2H14I ]. or, more recently, to study the relevance 
of thermal magnetic monopoles to the properties of the 
deconfined phase [15|-|20j. 

The detection of monopoles in lattice configurations 
is a highly non trivial problem. The procedure is well 
defined and gauge invariant in compact U(l) gauge the- 
ory [2l|: any excess over 2n of the Abelian phase of a 
plaquette is interpreted as existence of a Dirac string 
through the plaquette; a net magnetic charge exists in an 
elementary cube when a net number of Dirac flux lines 
crosses the plaquettes at its border. In this model the 
phase of the elementary Wilson loop is gauge invariant 
and therefore the procedure is unambiguous. In the case 
of a non-Abelian gauge theory, instead, one has first to 
fix a gauge, and then to apply the same procedure to the 
Abelian subgroup spanned by some diagonal component 
of the Lie algebra |3j. The result strongly depends on 
the choice of the gauge and, as a result, the existence 
of a monopole in a location of a given lattice configura- 
tion is a gauge dependent property, and this is of course 
physically unacceptable. 

In the soliton monopole of Refs. [22l l23j the Abelian 
subgroup which identifies the magnetic U(l) coupled to 
the magnetic charge coincides with the invariance sub- 
group of the vacuum expectation value of the Higgs field 



which produces the symmetry breaking. In QCD there 
is no Higgs field, but the magnetic (7(1) has to be a sub- 
group of the gauge group, due to the general property 
that the monopole component of any field configuration is 
intrinsically Abelian [24| . It was proposed in Ref. [3[ that 
any operator in the adjoint representation could be used 
as an effective Higgs field to identify the magnetic U(l) 
subgroup, the physics being for some reason independent 
of that choice. Each choice was called an Abelian pro- 
jection. In practice different Abelian projections proved 
to have different features, and there was a general con- 
sensus on the choice of the so-called "maximal Abelian 
gauge" [25| as the most convenient to expose Abelian 
dominance and monopole dominance (l2l . fl3t ] . More re- 
cent studies have shown that the differences between 
Abelian projections are less important than they look in 
simulations with low statistics [HHi], but the question 
of understanding the differences between Abelian projec- 
tions is still important. 

In a previous paper [23], which we shall quote as I in 
the following, these problems have been analyzed with 
the following results: 

1. The magnetic currents in any Abelian projection 
are proportional to the violation of the non-Abelian 
Bianchi identities, which is gauge covariant. They 
are the projection of that violation on the funda- 
mental weights which identify the Abelian projec- 
tions. 

2. Magnetic currents observed in lattice configura- 
tions depend on the Abelian projection. 

3. For each field configuration there exists a pre- 
ferred direction in color space, which is that of the 
(Abelian) magnetic monopole term in the multipole 
expansion [24]. That direction is the same that is 
identified by the maximal Abelian gauge. 

4. Only in the maximal Abelian gauge (modulo an ar- 
bitrary global transformation times a gauge trans- 
formation which is trivial at infinity) does the mag- 
netic charge obey the Dirac quantization condition. 
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In other gauges the charge measured by the flux of 
the magnetic field at large distances is generically 
smaller than the true monopole charge. 

5. Monopole condensation is a gauge-invariant fea- 
ture. 

Since the recipe of Ref. [2l| for counting monopoles 
is based on the Dirac quantization condition, all this is 
equivalent to state that the existence of a monopole is a 
gauge- invariant fact, but the detection recipe is Abelian 
projection dependent. 

All these findings were analytic and were analytically 
tested on the soliton solution of Refs. (22l. [23j. 

In the present paper we want to check them numeri- 
cally on Monte Carlo generated field configurations. As 
shown in Refs. |27ll2^|. proving this for SU{2) gauge the- 
ory in absence of fermions is equivalent to proving it for 
the generic case. We shall therefore simulate quenched 
SU{2). 

In Sec. II we show that the unitary gauge is noth- 
ing but the maximal Abelian gauge, as already stated in 
Ref. [271 . We then note that the "hedgehog" gauge of 
Refs. [22l[23| is nothing but the Landau gauge, and that 
in this gauge the flux at infinity of the magnetic field as 
defined by the 't Hooft tensor is vanishing. This explains 
why no monopole is found in the Landau gauge by use 
of the recipe of Ref. [U (see e.g. (H). 

We then introduce a class of gauge transformations 
continuously dependent on one parameter, connecting 
these two gauges, and we compute the expectation of 
the observed number of monopoles along this connec- 
tion, to be compared with the numerical data. The idea 
is that in the maximal Abelian gauge all the monopoles 
are detected, while, by changing the parameter of the 
gauge transformation, the effective magnetic charge is 
decreased and so is the number of observed monopoles, 
which is finally zero in the Landau gauge. 

Sec. Ill reports the numerical simulations and the re- 
sult of the check . 

In Sec. IV we discuss the results and we conclude. 



II. FROM MAXIMAL ABELIAN TO LANDAU 
GAUGE 

We consider the soliton solution of Refs. [13, [23j]. It 
was shown in the paper I that the result also applies to a 
generic configuration. With the notation of Ref. [3(| the 
gauge field of the soliton solution in the hedgehog gauge 
is 



Aa 
A 







A«=-(1-X(r))^- 
gr 



(1) 



where a is the color index, (0,n) the space time indexes, 
and K(r) a form factor vanishing for large r whose spe- 
cific form depends on the parameters of the Higgs sector. 



It is trivial to check that the solution in ([1} obeys the 
Landau gauge condition 



MS 



(2) 
The 't Hooft tensor in this gauge is by definition 



F^u = d^Al - d v A\ (3) 

and, by using the expressions in Eq. ([1]) and recalling 
that the Abelian fields are given by e% — Fqi and bi — 
TjCijkFjk, respectively, one gets for large r 



e = 0, 



„ T Z 

bi = -2 — 
gr 4 



(4) 



The magnetic charge Q which is detected by the recipe of 
Ref. [2l| is proportional to the flux of b across a spherical 
surface at r — >• 00, so one gets 



an (b-n) = - I an cos e = o (5) 
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No magnetic charge is thus expected in the Landau 
gauge. 

In the unitary gauge, instead, the field obeys the gauge 
condition which defines the maximal Abelian gauge, 
namely, 



(6) 



This can be easily checked by direct computation on the 
explicit solution given in the appendix of Ref. [30] . The 
exact field at all distances A^ in the unitary gauge, which 
is given in Eq. ([7]) below, exactly obeys Eq. (|6]): 



An = 



cos 6 — 1 




+ (1 - K ) sin(0 - 6) ) a 3 
sin 



01 



(cos(</>) + ia 3 sm((/>)) 



(7) 



Here = 9 °g+ c -i is a regulator of the singularity at 

9 = 7r and 0' its derivative with respect to 9; <p and 9 are 
the versors of the <fi and 9 directions. 

As shown in the paper I in this Abelian projection the 
Abelian magnetic charge is equal to two units, and can 
be detected by the recipe of Ref. 2l| since the Dirac 
condition is satisfied. The unitary transformation from 
the unitary (maximal Abelian) gauge to the hedgehog 
(Landau) gauge is known to be 



An = uA^W + Ud^u)W 



u 



(8) 

exp {-i<t>^) exp (-«©^) exp (i<f>^) 
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We shall operate a class of gauge transformations U(a) 
on the configurations in the maximal Abelian gauge, de- 
pending on one parameter a, such that for a = we 
stay in the maximal Abelian gauge [U (a) = 1] while for 
a = 1 we have U(a) — U of Eq. (|SJ) and we go to the 
Landau gauge. For each value of the parameter a, U (a) 
identifies an Abelian projection and the corresponding 
Abelian magnetic charge Q m {a) can be computed. We 
already know that Q m (0) = 2 and Q m (l) = 0. The 
general form of this transformation is 



U(a) =exp(-i7(0,0,a)y) x 

x exp (-ifi(9, (/), a) yj exp (ij(9, <p,a)^-) (9) 



with f3, 7 functions of the polar angles 9, <f> and of the 
parameter a, with boundary conditions 

7(0,0,1)=^ /3(M,0) = /3(0,<M) = 6 (10) 

We shall work out in detail the special case 

7(0, (f>,a) = 4> /3(9, <j), a) = a© (11) 

We are interested in the radial component of the 
Abelian magnetic field at large distances and along 03 
in color space. At large distances and 9 ^ tt Eq.© gives 
for the field 



■j* 1 7 cos 9 - 

A = -— (pa 3 — — 
Zgr smf 



1 



(12) 



Operating the transformation ^ and projecting on the 
third axis gives 



4" 
^3 



Tr 



0"3 



A a 



1 — cos0cos(a6 
2gr sin0 



and for the Abelian magnetic field 
1 



6(a) 



d 8 (sm9</)-A%) 



(13) 



(14) 



rsin9 

For the magnetic flux at infinity we have 

$(a) = r 2 y dfi f • 6(a) = - [l + cos(aTr)] (15) 

and for the ratio of the magnetic charge to that of the 
maximal Abelian projection we finally get 



i (a) 1 + cos(7ra) 



Qm(0) 



(16) 



This result is exact for the 't Hooft-Polyakov monopole 
and, by use of the same argument introduced in the paper 
I, it is also exact for a generic monopole in the continuum 
theory. 

On the lattice Eq. (Ti"6l) should give the ratio of the 
number of monopoles observed in the gauge correspond- 
ing to a generic value of a to that observed in the max- 
imal Abelian gauge. Of course on the lattice there are 



discretization errors which are not taken into account in 
the derivation of Eq. (TTB1) : a lattice configuration in the 
maximal Abelian gauge satisfies Eq. (j6|) only approxi- 
mately, the position of the monopole is not determined 
with arbitrary precision, the string direction is known 
with an angular error of vr/4, and the gauge field is con- 
stant along links. 

Because of these discretization errors we can not ex- 
pect Eq. (fl6)) to be exactly satisfied in lattice simula- 
tion. Nevertheless, we expect it to hold qualitatively: by 
increasing the parameter a from to 1 in a given con- 
figuration no new monopole should appear, but some of 
them should disappear, till the Landau gauge is reached 
at a = 1, where no monopole should be observed. 



III. NUMERICAL RESULTS 

As anticipated in Sec. I we shall test the above anal- 
ysis in quenched 577(2) since the results can trivially be 
exported to a generic gauge theory with and without 
quarks. 

To simulate SU (2) gauge theory we used a standard 
combination of heatbath [3l|, [32[ and overrelaxation [33[ 
updates. The maximal Abelian gauge fixing was achieved 
by an iterative combination of local maximization and 
overrelaxation steps (see e.g. the appendix of Ref. [34|]) 
and the algorithm was stopped when the average square 
modulus of the non diagonal part of the operator X(r), 
which has to be diagonal in the maximal Abelian gauge, 
was less than 10 -11 . Abelian links are then extracted 
as the third component of the gauge fixed field and 
monopoles are located by using the recipe of Ref. [2l[ . 

On this last point some care has to be used; for Eq. 
(ITBl to be true it is mandatory to integrate the Abelian 
magnetic flux on a surface at infinity, otherwise a residual 
dependence on the function K in Eq. (fT|), ([7]) is present. 
On the lattice monopoles are usually detected by inte- 
grating the flux on a cube whose faces are the elementary 
plaquettes. Is this sufficient to obtain the charges? It is 
possible to answer this question by using larger cubes, 
whose faces are n x n Wilson loops with n > 1: if the 
result is independent of n this means that elementary 
cubes are sufficient. The observed behavior is different 
in different gauges: while for the maximal Abelian gauge 
n = 1 is "large enough" (values obtained with n > 1 typ- 
ically differ less than 10% from the n = 1 ones), for the 
so-called local unitary gauges it is not [H, H(| ; we thus 
use just the elementary cubes in the gauge fixing, which 
is also what is done in the literature. We note, however, 
that this problem should worsen as we approach the con- 
tinuum limit: as j3 —> oo the physical scale grows and it 
is no more allowed to use n = 1. 

Once the monopoles are identified in the maximal 
Abelian gauge fixed configuration, the transformation (|9]) 
has to be applied, identifying the z axis with the direc- 
tion of the string which goes out of the monopole; we 
assume that the monopole is located in the center of the 
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cube. It is possible for the monopole to be located in- 
side a cube with more than one plaquette pierced by a 
string; if this happens it would be necessary to follow the 
strings and look for loops to recognize which is the one to 
be used as the z axis. For simplicity we restrict ourselves 
to the sample of monopoles with only one string, which 
should also be the great majority in the continuum limit. 
The gauge transformation © is thus applied to all the 
links of the elementary cube enclosing the monopole; the 
Abelian links are again extracted and the new magnetic 
charge inside the cube calculated. This charge obviously 
assumes only discrete values, so that to verify the relation 
(p~6|) we have to perform an average over the ensemble of 
the monopoles found in the configurations. 

Since in the entire procedure we have to operate only 
on the elementary cubes containing the monopoles, we 
can use a small lattice, of size 4 x 8 3 (the results were 
nevertheless checked for consistency also on a 4 x 16 3 
lattice). To compute mean values, 3 x 10 4 independent 
gauge configurations were generated at three different /3 
values: one below the deconfinement transition (f3 = 2.2), 
and two above the transition ((3 = 2.5 and {3 = 2.9). 




FIG. 1: Results and theoretical prediction (black line) for the 
ratio defined in (|16|) , 

The results of simulations are shown in Fig. ([1]) to- 
gether with the theoretical prediction Eq. (|T6| . We 
have to note that several sources of systematic error are 
present: 

1. we supposed the monopole to be in the center of 
the cube; 

2. the direction of the string is identified with an an- 
gular precision of only 7r/4; 

3. Eq.© is not exactly satisfied; 

4. the gauge field is defined only on the links and not 
pointwise. 



We have tested that the result is not much affected by 
the sources of error (1) and (2). The consequences of the 
lattice discretization (3) on the argument of the previ- 
ous section are hard to quantify. Nevertheless numerical 
results are in good qualitative agreement with the the- 
oretical expectations: for all nonvanishing a values the 
number of observed monopoles is strictly lower than the 
one in the maximal Abelian gauge and it approaches zero 
as a goes to 1 (Landau gauge). 

The deviations from the expected behavior and the 
residual /? dependence in Fig. Q] are probably due to lat- 
tice artifacts. Understanding in detail these artifacts 
goes beyond the scope of the present paper, whose main 
aim was to demonstrate that the gauge invariance of 
monopoles is not in contradiction with the gauge depen- 
dence of the number of monopoles observed in lattice 
configurations. 



IV. CONCLUSIONS 

In recent years evidence has been obtained that the 
existence of flux tubes is a feature not only of the max- 
imal Abelian gauge, but also other gauges, such as local 
unitary gauges, random Abelian projection gauge, and 
Landau gauge ( [3, HE HI, [S3] ) ■ These observations, to- 
gether with the ones obtained by looking for monopole 
condensation in Refs. [8H10|. strongly support the dual 
superconductor picture of gauge theories vacuum and the 
result of [27J that monopole condensation is a gauge- 
invariant phenomenon. 

In this work we have analyzed some of the conse- 
quences of the results derived in [27| for the numerical 
detection of monopoles in lattice configurations. The 
aim was to demonstrate the gauge invariance of the 
monopoles, an obvious prerequisite for any mechanism of 
confinement based on dual superconductivity. The use of 
non- Abelian Bianchi identities (see [27]]) allows us to de- 
fine monopoles in a gauge-invariant way and to explain 
why the observation of monopoles in lattice configura- 
tions does instead depend on the gauge. 

To test these results we considered a class of gauges in- 
terpolating between the maximal Abelian and the Lan- 
dau gauge and we calculated by a semiclassical model 
the expected variation in the number of the monopoles. 
Numerical results qualitatively agree with this theoreti- 
cal prediction. In particular we gave the first proof that 
in Landau gauge no monopoles can be detected by using 
the recipe of Ref . [2l[ , a long fact observed in numerical 
simulations but never understood. 
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